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Abstract 



A feasible mechanism of unitarization of amplitudes of deep inelastic scattering at small 
values of Bjorken x is the gluon fusion. However, its efficiency depends crucially on the vacuum 
color screening effect which accompanies the multiplication and the diffusion of BFKL gluons 
from small to large distances. From the fits to lattice data on field strength correlators 
the propagation length of perturbative gluons is Rc — 0.2 — 0.3 fermi. The probability 
to find a perturbative gluon with short propagation length at large distances is suppressed 
exponentially. It changes the pattern of (dif)fusion dramatically. The magnitude of the 
fusion effect appears to be controlled by the new dimensionless parameter ~ R'^/8B, with the 
diffraction cone slope B standing for the characteristic size of the interaction region. It should 
slowly oc l/lnQ^ decrease at large Q^. Smallness of the ratio R^/SB makes the non-linear 
effects rather weak even at lowest Bjorken x available at HERA. We report the results of 
our studies of the non-linear BFKL equation which has been generalized to incorporate the 
running coupling and the screening radius Rc as the infrared regulator. 
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1. Introduction. 

In processes of deep inelastic scattering (DIS) the density of BFKL [1] gluons, J-'{x,k'^), 
grows fast to smaller values of Bjorken x, T{x,k^) oc x~^, where, phenomenologically, A ^ 
0.3. The growth of T{x, k^) will have to slow down when the gluon densities become large 
enough that fusion processes gg ^ g become important. It was the original parton model idea 
of Refs. [2, 3] developed further within QCD in [4, 5]. The BFKL dynamics of saturation of 
the parton densities has been discussed first in [6, 7, 8], for the alternative form of the fusion 
correction see Eq.(AlO) of Ref. [9]. The literature abounds with suggestions of different 
versions of the non-linear evolution equation, see e.g. [10]. 

There is, however, at least one more mechanism to prevent generation of the high density 
gluon states. This is well known the vacuum color screening. The non-perturbative fluctu- 
ations in the QCD vacuum restrict the phase space for the perturbativc (real and virtual) 
gluons introducing a new scale: the correlation/propagation radius Rc of perturbative gluons. 
The perturbative gluons with short propagation length, Rc ~ 0.2 — 0.3 fermi, as it follows 
from the flts to lattice data on fleld strength correlators [11], do not walk to large distances, 
where they supposedly fuse together. The fusion probability decreases. We show that it 
is controlled by the new dimensionless parameter R^/SB, with the diffraction cone slope B 
standing for the characteristic size of the region populated with interacting gluons. 

The effects of finite Rc are consistently incorporated by the generalized color dipole (CD) 
BFKL equation (hereafter CD BFKL) [12, 13]. In presence of a new scale the saturation 
phenomenon acquires some new features and the goal of this communication is to present 
their quantitative analysis. 

2. CD BFKL and phenomenology of DIS. 

We sketch first the CD BFKL equation for qq dipole-nucleon cross section a{^,r), where 
^ = ln(xo/x) and r is the q'g-separation. The BFKL cross section cr(^,r) sums the Leading- 
Log(l/a;) multi-gluon production cross sections within the QCD perturbation theory (PT). 
Consequently, as a realistic boundary condition for the BFKL dynamics one can take the 
lowest PT order gg-nucleon cross section at some x — Xq. It is described by the Yukawa 
screened two-gluon exchange 
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(7(0, r) = (To(r) = 4Cf J _^ 2 \2 ^s{k'^)o:s{'i'^) 



x[l- Jo(A;r)][l-F2(k,-k)], (1) 

where — ^/Rc, Ois{K^) — 47r/^o ln(K^/AQ(^^) and — max{k^,C^ /r^}. The two-quark 
form factor of the nucleon can be related to the single-quark form factor 

F2ik,-k)^Fi{ue). (2) 

The latter is close to the charge form factor of the proton -Fi(g^) ~ Fp{q^) — 1/(1 q^/h?Y, 
where = 0.71 GeV^ and in Eq.(2) u = 2NJ{N^ - 1) for the color group SU{N^) [14]. 

The small- a; evolution correction to a{^, r) for the perturbative 3-parton state qqg is as 
follows [12, 13] 

9^a(e,r) = Jd'p, |V(Pi)-V'(P2)l' 

X [(^3(C,r,Pi,p2) -t^(C,^^)], (3) 

where the 3-parton (gg^f-nucleon) cross section is 

C 

(73(e, r, pi, P2) = ^ [a(e, Pi) + cT(e, P2) - r)] + a{C, r), (4) 

where Ca — Nc and Cp — {N^ — l)/2Nc. Denoted by Pi^ 2 the q-g and q-g separations 
in the two-dimensional impact parameter plane for dipoles generated by the q-q color dipole 
source. The radial light cone wave function V'(p) of the dipole with the vacuum screening of 
infrared gluons is [12, 13] 



m = ^^^^^^-^K,{p/R^l (5) 
where Ki,{x) is the modified Bessel function. The one-loop QCD couphng 

as{Ri) = 47r//3o HC^/KIcdRD (6) 

is taken at the shortest relevant distance Ri = min{r, pi}. In the numerical analysis C = 1.5, 
^QCD — 0.3 GeV, Po — {llNc — 2Nf)/3 and infrared freezing as{r > rf) — af — 0.8 has 
been imposed (for more discussion see [15]). The scahng BFKL equation [1] is obtained from 

Eq. (3) at r, pi^2 ^ Rc in the approximation as = const - the dipole picture suggested in [16]. 



3. Perturbative and non-perturbative. 

The perturbative gluons are confined and do not propagate to large distances. Available 
fits [11] to the lattice QCD data suggest Yukawa screening of perturbative color fields with 
propagation/screening radius Rc ~ 0.2 — 0.3 fm. The value Rc — 0.275 fm has been used since 
1994 in the very successful color dipole phenomenology of small-x DIS [17, 18, 19, 20, 21]. 
Because the propagation radius is short compared to the typical range of strong interactions 
the dipole cross section obtained as a solution of the CD BFKL equation (3) would miss the 
interaction strength for large color dipoles. In [17, 18] this missing strength was modeled by 
the ^-independent dipole cross section and it has been assumed that the perturbative, (t(^, r), 
and non-perturbative, (7„pt(r), cross sections are additive, 

(^totii, r) = a{i, r) + anpt{r). (7) 

The principal point about the non-perturbative component of (Ttot{,i,r) is that it must not 
be subjected to the perturbative BFKL evolution. Thus, the arguments about the rise of 
(7(^, r) due to the hard-to-soft diffusion do not apply to (7„pt(r). We reiterate, finite Rc means 
that gluons with the wave length \ Rc are beyond the realm of perturbative QCD. A 
quite common application of purely perturbative non-linear equations [6, 7] to the analysis 
of DIS data without proper separation of perturbative and non-perturbative contributions is 
completely unwarranted. 

Specific form of (Jnpt{T) motivated by the QCD string picture and used in the present paper 
is as follows: 

(Tnpt{r)^aal{r)r''/{r + d). (8) 

Here d = 0.5 fm is close to the radius of freezing of the running QCD coupling r/ and a = 5. 
fm. 

Our choice Rc — 0.26 fm leads to a very good description of the data [22, 23, 24, 25, 26] 
on the proton structure function F2{x,Q^) at small x shown in Fig.l. Shown separately are 
the nonperturbative contribution (8) and the contribution from DIS off valence quarks [27]. 
The effects of quark masses important at low Q'^ are taken into account [28]. The linear 
CD BFKL description of F2{x, Q"^) (dashed hue) is perfect at moderate and high Q"^ where 
it is indistinguishable from the solid line representing the non-linear CD BFKL results (see 
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below) . Two lines diverge at low Q"^ where the account of the non- linear effects improves the 
agreement with data. 

Recently a global analysis of HERA DIS data has been reported [29]. In [29] a purely 
perturbative non-linear equation is solved with some phenomenological initial conditions. A 
very soft infrared regularization with the infrared cutoff ~ ^c^cd allows non-perturbatively 
large dipoles to be governed by the perturbative QCD dynamics. The non-perturbative 
component of solution evolves perturbatively to smaller x. Good agreement with data was 
found. 

4. CD BFKL and the partial-wave amplitudes. 

Following [30, 31] we rewrite the Eq.(3) in terms of the g'^-nucleon partial- wave amplitudes 
(profile functions) r(^,r, b) — 1 — 5'(^,r, b) related to the scattering matrix 5'(^,r, b). We 
introduce the impact parameter b defined with respect to the center of the q-q dipole. In 
the qqg state, the qg and qg dipoles have the impact parameter b -|- P2,i/2. In the large- A^c 
approximation 173 in Eq. (4) reduces to (73 = cr(C,Pi) + <^(C,P2)- what corresponds to the 
factorization of the 3-parton {qqg) scattering matrix, 

SsiC, r, Pi, P2) = S{C, pi, b + ^P2)S{C, p2, b + ^p,). (9) 

Then, the renormalization of the g'^-nucleon scattering matrix, S{^, r, b), for the perturbative 
3-parton state qqg is as follows 

(10) 



X 



S{t pi, b + ^P2)^(e, P2, b + ^p,) - 5(e, r, b) 



For the early discussion of Eq. (10) see [6, 7]. The substitution 5'(^, r, b) — 1 — r(,^,r,b) 
results in 

9^r(e,r,b)= Jd'p,mp,)-^l;(p,)f 

X r(e, pi, h + ^p,) + r(e, P2, b + ^pi) - r(e, r, b) 

-r(e, P2, b + ^Pi)r(e, pi, b + ^P2). (11) 
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We identify the corresponding partial waves using the conventional impact parameter repre- 
sentation for the elastic dipole-nucleon amplitude 



;i2) 



For the predominantly imaginary /(i^,r, k) = ia{^,r) exp{—Bk^ /2) the profile function is 



r(e,r,b) = /^^;;\ exp 

47rS(^,r) 

anda(e,r) = 2/rf2br(e,r,b). 

Integrating over b Eq. (11) yields [32] 



2S(e,r), 



(13) 



d^aitr)^ Jd'p, |V'(Pi)-V'(P2)r 



47r(5i + B2) 



exp 



8(5i + B2 



(14) 



where Bi — B{^, pi). The above definition of the scattering profile function, Eq. (13), removes 
uncertainties with the radius R of the area within which interacting gluons are expected to 
be distributed (the parameter Sx_ = t^R^ appearing in Eq. (25)). In different analyses of the 
non-linear effects its value varies from the realistic R? — 16 GeV~^ [33] down to the intriguing 
small i?^ = 3.1 GeV~^ [34]. Besides, the radius R is usually assumed to be independent of x. 
In our approach the area populated with interacting gluons is proportional to the diffraction 
cone slope 5(^, r). 

5. The diffraction cone slope. 

The diffraction slope for the forward cone in the dipole-nucleon scattering [30] was pre- 
sented in [31] in a very symmetric form 



B{Lr) = \{h'')^\T' + \Rl + 2a'^i. 



(15) 



The latter provides the beam, target and exchange decomposition of B: is the purely 
geometrical term for the color dipole of the size r, Rn represents the gluon-probed radius of 
the proton, the dynamical component of B is given by the last term in Eq. (15) where a'jp is 
the Pomeron trajectory slope evaluated first in [30] (see also [31]). The order of magnitude 
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estimate [31] 

^'^-^j as{r)Rfr'Kf{r/R,) ^ ^asiRc)Rl , (16) 

clearly shows the connection between the dimensionful a'jp and the non-perturbative infrared 
parameter Rc. The increase of B with growing collision energy is known as the phenomenon 
of shrinkage of the diffraction cone. 

We determine a'jp as the ^ — >■ oo limit of the local Regge slope a'^ff{^,r) — d^B{^,r)/2 
[31]. At ^ — )■ oo, a'^ff{^,r) tends to a r-independent a'jp = 0.064 GeV~^. The onset of the 
limiting value a'jp is very slow and correlates nicely with the very slow onset of the BFKL 
asymptotics of a{^, r) [12]. An interesting finding of Ref. [31] is a large sub-asymptotic value 
of the effective Regge slope a'^jj{^, r), which is by the factor ~ (2 — 3) larger than a'jp. 

In Eq. (15) the gluon-probed radius of the proton is a phenomenological parameter to be 
determined from the experiment. The analysis of Ref. [35] gives i?^ w 12GeV~^. 

6. Non-linear CD BFKL: small dipoles, r < R^. 

The term quadratic in a in Eq. (14), models the process of the gluon fusion. The efficiency 
of this "fuser" differs substantially for r ^ Rc and for r ^ Rc- Consider first the ordering of 
dipole sizes 

< < Rl (17) 
corresponding to the Double Leading Log Approximation (DLLA) [36]. Eq.(14) reduces to 



TT 



dp 

X ' 



(18) 



SttB 

First notice that the function 

g{C,p)^p-MC,p) (19) 
is essentially flat in p and the second term in the rhs of Eq. (18) is dominated by p ~ R^ 



1 „^2., ^: 

p' 

with B = B{^,Rc). Thus, the new dimensionless parameter 



1 r^c dn^ R^ 



enters the game. Its geometrical meaning is quite clear. Remind that the unitarity requires 
(see Eq. (13)) 

(7(e,p)<87rS. (22) 

Smallness of k,c makes the non-linear effects rather weak at HERA even at lowest available 
Bjorken x (see Fig. 1). Comparison of the linear and quadratic terms in the right hand side of 
Eq. (18) shows that the relative strength of non-linear effects decreases to smaller ~ Q~'^ 
logarithmically 

K=^^o,K,ln-\Q'Rl), (23) 
lin. 

Therefore, we are deahng with the scahng rather than the higher twist, effect. 
7. Saturation scale and observables. 

The parameter k in Eq. (23) should not be confused with another parameter frequently 
used to quantify the strength of the non-linear effects. It decreases with growing much 
faster than k in Eq. (23). Namely, 



K ~ — . (24) 



1 

The estimate (24) comes from equating the linear and non-linear terms in rhs of the equation 

[4, 5] 

d^d.Gi^, rj) = cG(e, v) - ^G(e, vf, (25) 

where r] = \n{Q^ / Aq^^), c = asNc/n, a = a'g7i/S±, G{^,ri) is the integrated gluon den- 
sity and Eq.(25) comes from Eq.(18) as for small dipoles g{^,p) ~ j^as{p)G{^, p). The 
corresponding value of Q'^ denoted by 

Ql^aG{x,Ql)/c (26) 

is called the saturation scale. The non-linear saturation effects are assumed to be substantial 
for all Q <:,Qs (see e.g. [37]). Obviously, Eq. (23) asserts something different. The point is that 
Eqs. (23) and (24), describe the Q^-dependence of strength of the non-hnear effects for two 
very different quantities: the integrated gluon density G{^, Vj) and the differential gluon density 
•^(C) v) — 9r)G{^, rj), respectively. The gluon density G(^, 77) is a directly measurable quantity. 
For example, the longitudinal DIS structure function is Fl(x, Q^) ~ as{Q'^)G{x, Q"^) [38]. On 
the contrary, the differential gluon density J-'({, Vj) is related to the observable quantities like 



F2{x, Q"^) rather indirectly, by means of the well known transformations leaving a weak trace 
of Eq. (24) in F2{x,Q^). 

A possibihty to test Eqs. (24) and (26) provides the coherent diffractive dijet production 
in pion-nucleon and pion- nucleus collisions [39]. Both helicity amplitudes of the process 
are directly proportional to J-'{x,k'^). The same proportionality of diffractive amplitudes to 
T{x,k'^) was found for real photoproduction with pointlike jqq vertex in [40]. Therefore, 
there is no real clash between Eqs. (23) and (24). The sharp Q^-dependence of the nonlinear 
term in Eq. (25) does not imply vanishing non-linear effects in G{x, Q"^) for Q"^ ^ Q^. 

8. Non-linear CD BFKL: large dipoles, r ^ Re- 

The interplay of the color screening and gluon fusion effects at large r ^ Rc, where the 
non-linear effects are expected to be most pronounced, requires special investigation. In high- 
energy scattering of large quark-antiquark dipoles, r ^ Rc, a sort of the additive quark model 
is recovered: the (anti)quark of the dipole r develops its own perturbative gluonic cloud and 
the pattern of diffusion changes dramatically. Indeed, in this region the term proportional to 
Ki{pi/ Rc)Ki{p2/Rc) in the kernel of Eq. (3) is exponentially small, what is related to the 
exponential decay of the correlation function (the propagator) of perturbative gluons. Then, 
at large r the kernel will be dominated by the contributions from pi <Rc^P2 — t and from 
P2 ^Rc^ Pi — f- It does not depend on r and for large Nc the equation for the dipole cross 
section reads 

asC, 



9MC,r) = ^ / d'p,RfK'MR,) 



{(7(e,pi)+a(e,P2)-(7(e,0 
o-(C,Pi)o-(C,P2) 



47r(Si + B2) 



exp 



^2 



S{Bi + B2) 



(27) 



where — B{^, pi). For a qualitative understanding of the role of color screening in the non- 
linear dynamics of large dipoles we reduce Eq. (27) to the differential equation. First notice 
that the dipole cross section cr(^, r) as a function of r varies slowly in the region r ^ Rc, while 
the function Ki{y) vanishes exponentially at y » 1 and Ki{y) f=:i 1/y for y <^ 1. Therefore, 
Eq. (27) can be cast in the following form 

c'%a{i, r) = a{i, R^) + Rld^i, r) - 
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-a{i,Rc)a{i,r)/8nB, (28) 

where c — asCp/'K and for simplicity B — B{0,Rc). The solution of Eq. (28) with the 
boundary condition cr(0,r) = cro(r^), where cro(r^) comes from Eq. (1), is 

, _ ao(r' + <fi;)+»K) 

l + v(i)/8nB ' ^ ' 

where 

i;(^) = e"^ ao{Rl + zRl)e-'dz. (30) 
io 

Prom Eq. (1) it follows that at r < Z = u^Yn.{RJ^/2, ^/u/K} 

Mr') - i^r^„.(r/VI) In (31) 

where ^4 fa 10 comes from properties of the Bessel function Jo{y) in Eq. (1) [41]. For large 
dipoles cro(r^) saturates at ~ AP, 

(7o(r^) ^ 47rCFRlafas{Rc)h{a), (32) 

where a/ = 0.8 (see Eq.(6)) and the interplay of two scales, Rc and A~^, in Eq. (1) results in 
h{a) — 1 — (a^ — 1 — 2a In a)/(a — 1)^ with a — it/i?^A^. This kind of saturation is due to the 
finite propagation radius of pcrturbative gluons. 

With growing ^ the dipole cross section cr(^, r) increases approaching the unitarity bound, 
(7 = SttB. To quantify the strength of the non- linear effects we introduce the parameter 

K — Sa/a, (33) 

where Sa — a — Uni and a represents the solution of the linear CD BFKL Eq. (3), while 
Gni stands for the solution of the non-hnear CD BFKL Eq.(14). Therefore, our k gives the 
strength of the non-linear effects with the non-perturbative corrections switched off 

K = oa/a 



v{^) + SnB 



~ 47rAve^a5(i?c/VI)(e^« - 1). (34) 
Po 

The magnitude of non-linear effects is controlled, hke in the case of small dipoles, by the ratio 

Rl/B (we assumed Rl -C u/A^). 
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Numerical solution of Eqs.(3,14) gives the r-dependence ol k = 5a / a shown in Fig. 2 for 
several values of a; and for two correlation radii Rc = 0.26 fm and Rc = 0.52 fm. For 5a /a <^ 1 
the law 5a/a (X Rl/B holds true. At large r ^Rc the toy- model solution, Eq. (29), (dashed 
lines) correctly reproduces the ^-dependence of k. At small r <^ Rc the ratio 5a /a decreases 
slowly as it is prescribed by Eq. (23). In Fig. 2 also shown is the evolution of the unitarity 
ratio aniit r)/(47r(Bi + B2)) with Bi ^ B{^, R^), B2 = -B(^, r) and denoted by a/8nB. High 
sensitivity of ani{C, t^) to Rc is not surprising in view of the toy-model solution (29). 

9. Summary. 

To summarize, the purpose of the present paper has been an exploration of the phe- 
nomenology of saturation in diffractive scattering which emerges from the BFKL dynamics 
with finite correlation length of perturbative gluons, Rc- The non-linear effects are shown to 
be dominated by the large size qq — g fluctuations of the probe (virtual gauge boson). They 
should very slowly, oc 1/ In Q^, decrease at large Q^. The magnitude of the non-linear effects is 
controlled by the dimensionless parameter Kc — R^/SB. The area populated with interacting 
gluons is proportional to the diffraction cone slope B. Smallness of Kc makes the non-linear 
effects rather weak even at lowest Bjorken x available at HERA. The linear BFKL with the 
running coupling and the infrared regulator Rc — 0.26 fermi gives very good description of 
the proton structure function F2{x, Q^) in a wide range of x and Q^. 
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Figure 1: The CD BFKL description of the experimental data on F2{x, Q^). Black triangles 
and circles are ZEUS data [22, 23], open triangles and circles show HI data [24, 25] and 
open squares refer to E665 results [26] . Dashed hues represent the linear CD BFKL structure 
function F2. Shown by sohd hues are the non-hnear CD BFKL structure functions F2. At 
high the non-linear effects vanish and both dashed and solid lines are indistinguishable. 
The valence and non-perturbative corrections are included into both the CD BFKL and the 
non-linear CD BFKL description of F2. The contribution to F2 from DIS off valence quarks 
[27] is shown separately by dash-dotted lines. Shown by dotted hues are the non-perturbative 
contributions to F2. 
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Figure 2: The dipole size dependence of the non-hnear correction k, — Sa / a to the hnear CD 
BFKL dipole cross section a for two correlation radii Rc and for ^ = 6, 8.5, 11, 13, 15.5, 20. 
Dashed lines correspond to the approximation v{^) a; aoiRDie"^^ — 1) in Eqs. (29,30). Shown 
separately is the "unitarity ratio" a/SnB (see text) for two values of Rc and for the same set 

ofC- 



16 



